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ABSTRACT 

We apply a theory of wave propagation through a turbulent medium to the scattering of radio 
waves in pulsar magnetospheres. We find that under conditions of strong density modulation the 
effects of magnetospheric scintillations in diffractive and refractive regimes may be observable. 
The most distinctive feature of the magnetospheric scintillations is their independence on 
frequency. Results based on diffractive scattering due to small scale inhomogeneities give a 
scattering angle that may be as large as 0.1 radians, and a typical decorrelation time of 10~ 8 
seconds. Refractive scattering due to large scale inhomogeneities is also possible, with a typical 
angle of 10 -3 radians and a correlation time of the order of 10~ 4 seconds. Temporal variation 
in the plasma density may also result in a delay time of the order of 10 -4 seconds. The different 
scaling of the above quantities with frequency may allow one to distinguish the effects of 
propagation through a pulsar magnetosphere from the interstellar medium. In particular, we 
expect that the magnetospheric scintillations are relatively more important for nearby pulsars 
when observed at high frequencies. 

1. Introduction 

A number of observational results may possibly be attributed to scattering processes inside pulsar 
magnetospheres. After most of the present work was completed, the very interesting and convincing 
results of Sallmen et al. (1999) were published. In this work, the authors found frequency independent 
(measured at 1.4 and 0.6 GHz) jitter in the arrival time of giant pulses from the Crab (of order 100/xsec). 
In addition, the frequency independence of the spread and the multiplicity of the pulse components with 
large variations in the pulse broadening times strongly suggests that multiple components of the giant 
pulses are due to refractive scattering inside the pulsar magnetosphere. Similarly, Hankins & Moffett (1998) 
found that broadening times for a single giant pulse from the Crab pulsar scale more slowly with frequency 
than A 4 (the value predicted if the scattering is due solely to the unmagnctized electron-ion interstellar 
plasma). Gwinn et al. (1997) determined, using interstellar scintillation, that the size of the Vela pulsar's 
radio emission region is about 500 km. This, being of order 1/10 the light cylinder radius, is considerably 
larger than conventional estimates. Similar results were obtained by Smirnova et al. (1996) and Wolszcan 
& Cordes (1987) (but see also Cordes, Weisberg k Boriakoff 1983). Gwinn et al. (1999) found that 
pulsar 0437-471 shows frequency independent scintillations with bandwidth 4 MHz at several observing 
frequencies below 1 GHz. Kramer et al. (1997) found that intensity fluctuations at very high frequencies 
(30 GHz) are considerably enhanced if extrapolated from from observations at lower frequencies. Kramer 
et al. (1999) found frequency-dependent changes in one of the pulse components of the msec pulsar PSR 
J1022+1001 these take place on a characteristic frequency scale of 8 MHz at the frequency 1410 MHz. With 
a dispersion measure of only 10 pc/cm 3 , this pulsar is probably too close to show diffractive interstellar 
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scintillations with such bandwidth. Citing the authors of the above reference "a propagation effect in the 
pulsar magnetosphere might be still the most probable explanation for the observed phenomena." 

Simple models for the interstellar medium fail to fully explain the scattering properties of nearby 
pulsars. This scattering is enhanced, relative to that extrapolated from more distant pulsars (Sutton 1971), 
and shows stronger refractive effects (Gupta et al. 1994, Bhat et al. 1999, Rickett et al. 1999). The 
proposition by Hajivasiliou (1992) and Bhat et al. (1998) that the excess scattering of nearby pulsars lies 
at the surface the Local Bubble was criticized by Britton et al. (1998) who found strong upper limits 
on the angular broadening of 5 nearby pulsars, and found that these sizes were inconsistent with both 
a uniform medium and scattering at the surface of the Local Bubble. Through the examination of the 
angular and temporal broadening of nearby pulsars, Britton et al. (1998) also concluded that the scattering 
material of nearby pulsars lies near the pulsars, and is moving in the same general direction as the pulsars. 
Rickett et al. (1999) found that scattering material for PSR 0809+74 probably lies close to the pulsar, 
using a technique based on the time scale of scintillation, although it could lie much closer if its velocity 
were aligned with that of the pulsar. Scattering in pulsar magnetospheres could contribute to interstellar 
scattering and eliminate the necessity of postulating a different, nearby, scattering medium. Taken together, 
these results compelled us to investigate the scattering of radio waves inside the pulsar magnetosphere. 

A theory of interstellar scattering is well developed and has been successfully applied to explain 
various effects observed in pulsars (Blandford & Narayan 1985, Lee and Jokipii 1975a, b,c, Rickett 1977). In 
this work we will study the effects of scattering and diffraction inside the pulsar magnetosphere. Strong 
magnetic fields present in pulsar magnetospheres and the unusual electrodynamics of the one-dimensional 
electron-positron plasma both change the familiar effects of scattering and refraction in plasma. The 
unusual features of scattering in such plasma may allow separation from the interstellar scattering and will 
serve as a tool to probe the structure of the magnetosphere itself. 

A large variety of wave-plasma interactions, including scattering and diffraction, may be described in 
terms of the variation of the refractive index of a medium. A superstrong magnetic field plays a dominant 
role in defining the electron-photon interactions at small energies. Typically, the frequency of the observed 
radio waves is much less that the cyclotron frequency uj <C u>b- In this limit, the strong magnetic field 
suppresses the wave-plasma interaction by (uj/usb) 2 (the Thompson cross-section is smaller by this ratio 
compared to an unmagnetized medium). For such frequencies the refractive index of strongly magnetized 
pair plasma for escaping electromagnetic modes is approximately 

n 2 -i«4, a) 

which is fundamentally different from the case of unmagnetized electron- ion plasma (where n 2 — 1 = — -f ) . 
In Eq. (Q) ujp is the plasma frequency, ujb is the nonrelativistic cyclotron frequency. Equation (Q) assumes 
that the plasma is stationary in the pulsar frame. On the open field lines, the pair plasma is moving with a 
larger streaming Lorentz factor 7. Then, the refractive index for waves propagating at comparatively large 
angles to the magnetic field 9 » 1 /j p is (see appendix |X|) 

n 2 -l«^*. (2) 

Thus, the large Lorentz factor of the moving plasma effectively enhances the wave-plasma interaction on 
the open field lines. 

Parameter S is the key to the scattering and diffraction effects in the pulsar magnetosphere. An 
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important fact is that 5 does not strongly depend on our assumptions about the density and the streaming 
Lorentz factors of the plasma. To see this, we normalize the density of the secondary plasma to the density 
of the beam, lo^ = XmoSqj — 2AmWbO, where Am is a multiplicity factor, and use the fact that the energy 
in the plasma is approximately equal to the energy in the beam, 7 p Am ~ lb- We find then 

6=™ (3) 
ujb 

Theoretical estimates (e.g. Arons 1983) of the beam's Lorentz factor give 76 w 10 6 ~ 7 . This value is 
determined by the condition of vacuum breakdown due to pair production, and depends mostly on the 
structure of the magnetic field (radii of curvature) in the acceleration zone. In the first approximation, 75 
may be considered a constant for all the pulsars, regardless of their spin period or surface magnetic field. 
Then for a conservative estimate 76 w 5 x 10 5 , parameter 5 near the light cylinder (where scattering is most 
important - see below) is 

This is a comparatively large value. For example, a fluctuation of the order of unity in the relative density 
inside the magnetosphere will induce time delays of the order of SP w 30/i sec; this is a small value if 
compared with the total interstellar time delay for typical pulsars, but it is comparable to the relative 
dispersive interstellar time delay between the frequencies separated by 1 GHz. 

Parameter 5, which determines refractive properties of the medium, is negligible deep inside the 
pulsar magnetosphere, but increases with the distance from the neutron star as cx r 3 . Thus, the strongest 
nonresonant wave-plasma interactions occur in the outer regions of pulsar magnctospheres (near the light 
cylinder). This allows for a considerable simplification when considering scattering and diffraction effects 
since one can adopt a "thin screen" approximation. We assume that emission is generated deep in the 
pulsar magnetosphere and then scattered in a thin screen located near the light cylinder with a typical 
thickness D w O.IRlc- 

The physical picture of the disturbances that we have in mind consists of small scale and large 
scale turbulence. The former is excited possibly by microscopic plasma turbulence or nonlinear plasma 
interactions (like Langmiur collapse or modulation instability), while the latter is mostly due to the 
temporal and spatial modulation of the outflowing pair plasma at the moment of its creation. We suppose 
that the typical sizes of the small scale inhomogeneities should be comparable to tens of Debye radii 
fl min ~ 10 x r D = l§Lo p /v T fts 3 x 10 2 cm (v T is a typical thermal velocity, r D is a Debye radius); while 
typical sizes of the large scale inhomogeneities are of the order of a tenth of the light cylinder radius 
a max ~ O.IRlc — 5 x 10 7 cm. As we will show later, the small scale inhomogeneities will contribute to 
diffractive scattering, while larger scale inhomogeneities will contribute to refractive scattering. 

Since a detailed spectrum of the density perturbations inside the pulsar magnetosphere is presently 
not known, we will also make estimates for the power law spectrum extending from smallest scales a m i n to 
largest scales a max . Relating these two qualitatively different effects by a simple power law may be less 
justified here than for the case of interplanetary and interstellar media. However, it allows for a consistent 
treatment of refractive and diffractive effects and, permits the use of well developed methods. We expect 
the results to depend qualitatively on either the inner or the outer scales and only weakly on the particular 
choice of the power law index a. 

Another difference from the case of interstellar scattering worth mentioning here is that while plasma 
in the pulsar magnetospheres is one dimensional, we still expect the small and large scale turbulence to be 
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three dimensional. At large scales this is an obvious consequence of the temporal and spatial modulation of 
the outflowing pair plasma; at small scales, the excitation of transverse turbulence (modulating densities of 
plasma across the magnetic field) is thought to be as effective as excitation of longitudinal turbulence (e.g., 
Weatherall) . 

Contradictory estimates exist regarding the development of small scale turbulence in the pulsar plasma. 
The primary source of this turbulence is usually associated with electrostatic instabilities which may not 



have enough time to develop ( Melrose 1995 ). On the other hand, the existence of the large scale density 
fluctuations is certain. If the small scale turbulence, which is responsible for the diffractive scattering (see 
below) , does not develop, then our estimates for refractive-type effects will still be valid []. 

There is an important unresolved problem which may affects the results of this work. Electromagnetic 
waves should be absorbed in the outer parts of the magnetosphere at the cyclotron resonance: 

uj S3 ujb/Ip- (5) 

This resonance occurs in the outer regions at 

^3xlO«c m m- 1/a .rir I/a feV 1/I (6) 



vio 9 / vioo/ V 1012 , 

Nonetheless, radiation avoids being absorbed. Q 

In addition, dispersion relation Eq. (|2|) is valid only for lo <C ujs/jp- In order to avoid the absorption 
at the cyclotron resonance, we limit our considerations to fast pulsars with periods P 

For numerical estimates, we will use pulsars with a period P — 0.1 sec, surface magnetic field 
£?ns = 10 12 G, plasma frequency on the open field lines cu 2 = 2 Am £Iujb, multiplicity factor Am = 5 x 10 3 
and a streaming Lorentz factor 7 P = 100. Then at the light cylinder radius Rlc = 5 x 10 s cm, the 
parameter for both open and closed field lines is S S3 3 x 10~ 4 . 



2. Order of magnitude estimates 

2.1. Diffractive and refractive scattering in the case of inhomogeneities of two scales 

Consider a screen of thickness D, with irregularities of typical size a. A variation in the refractive index 
An extending over a length a induces a change in the phase of a wave by 8<j) ss waAn/c. Using dispersion 
relation Eq. (^|) the variation in the refractive index may be written as 

An S3 4nr e c 2 ^ p AN e /uj% = SAn e (8) 



1 In the absence of diffractive scattering, refractive-type scattering may result in a formation of caustics. 

2 The fact that we do see pulsar radio emission implies that the dipole approximation breaks down in the outer magnetosphere. 
This may be interpreted in two ways: (1) there is a sudden decrease in the density of plasma before the cyclotron resonance, 
e.g., due to the sweepback near the leading last open field line, (2) the field falls off slower than 1/r 3 . 




Fig. 1. — The emission is generated deep in the pulsar magnetosphere and then scattered in a thin screen 
located near the light cylinder with a typical thickness D « O.IRlc- 
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where r e = e 2 /mc 2 is the classical radius of an electron, and we have normalized density perturbations as 
An e = AN e /N e . The phase shift relative to a vacuum-propagating wave after propagating the distance a, is 

S(f> i=a 2ttS An e ^ oc A -1 . (9) 
A 

The inverse wavelength dependence of the phase shift highlights the important point: scattering becomes 
stronger at high frequencies. This is opposite to the scattering trend by ISM. 

A ray passing through the whole screen encounters D/a randomly distributed irregularities; therefore, 
the difference between the phases of rays separated laterally by more than a is 

m i/2 /m 1/2 .. a 



a 



A(j)^2n[ — ) 5<f>=[ — ) 8An e -. (10) 



The form of the scattered image will depend on whether the scattering is weak (when S(f> <C n) or 
strong (when 5<p^> -n). The angular spectrum of waves has two limiting cases: 

9 S = -, A</> < 7T 

a 

e s = ^-, Ac^yyn (11) 
Z7r a 

In the case of weak scattering, we see the initial wave and a thin halo around the mostly unscattered 
wave. In this limit, the transverse dimension of the correlation function, p c — \/6 Cl is equal to the size 
of inhomogeneities. In the case of strong scattering, most of the energy of the initial wave has been 
transformed into a scattered component. The transverse dimension of the correlation function becomes 
smaller: p c w a/A<p. 

In the case of scattering inside the pulsar magnetosphere, the typical phase delay A(j) is larger than tt 
(see below) - this implies that scattering occurs in the strong regime. The rays are typically scattered by 
the angle 

^«^=f-) 1/ %An e . (12) 
2-k a \ a J 

The two types of inhomogeneities that should be present inside the pulsar magnetosphere will produce 
qualitatively different effects: small scale inhomogeneities will produce diffractive scattering, while large 
scale inhomogeneities will produce refractive scattering. To see this, we recall (e.g. Prokhorov et al. 1975) 
that strong scattering may occur in refractive, diffractive or mixed regimes. Separation into refractive 
and diffractive effects is based on the observation that the geometric optics limit becomes unsuitable for 
distances larger than 

L CI =^. (13) 

For larger distances, the diffractive pattern due to the refractive inhomogeneities becomes important 
regardless of the small size of the diffraction angle. Thus we can identify two scattering regimes 
(i) diffractive, where the inhomogeneities are seen mostly in the Fraunhofer zone and (ii) refractive, 
where the inhomogeneities produce fluctuations in the index of refraction and can be described by the 
geometrical optics approximation. The two limits are separated by the size of the first Fresnel zone, 
ff = \/Ai?LC ~ 10 5 cm. Inhomogeneities with a size less than the size of the first Fresnel zone contribute 
to diffractive scintillations, while those with a size larger than the first Fresnel zone contribute to refractive 
scintillations. 
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For diffractive scattering, Eq. ( |12| ) is an estimate of the the apparent angular size of the source seen 
through the screen. Numerically for D ss Q.IRlc, a — a mm and under the most optimistic assumption of 
relative density fluctuations of the order of unity An e sa 1, we find: 

D » O.lrad. (14) 

This implies that large angle scattering is possible in the outer regions of pulsar magnctosphcrcs. The 
observed profiles are then the convolution of the "initial" window function (determined by the emission 
conditions at lower radii) with diffractive scattering. 

We would like to reiterate here the the estimate of the scattering angle in Eq. (|l4|) assumes that all the 
power in the density fluctuations is concentrated at one diffractive scale and that the turbulence is strong 
< AiV e >« N e . Given the total power in fluctuation (e.g. < An e > 2 w 1), the assumption of a power 
law spectrum of density fluctuations removes some power from small scales and reduces the diffractive 
scattering angle (see Sec. 2.2). On the other hand, parameter S increases with period of the pulsar - in the 
longer period pulsars it will be larger than we assumed. 

The diffractive image will be focused and defocused by the refractive fluctuations with a scale 
« OdRlc ~ o-max- A typical refractive scattering angle is 

6> R «<S = 4x 1(T 3 . (15) 



Refractive effects will induce "jitter" in the arrival times of the pulses and a temporal correlation in 
the intensities with a typical scale tr = 6*rP = 4x 10 _4 sec. Both of these effects will be independent of 
frequency, and increasing with the period of the pulsar. 

As we have already mentioned, a necessary requirement to allow the separation of scintillations into 
diffractive and refractive branches is that the scattering should be strong (the total phase shift (Eq. (|l0|)) is 
much larger than tt). Using our simple estimates (Eq. (|l0|)) for a medium with given size inhomogeneities, 
we find that this condition is satisfied for 

A 2 

a win >^=2x 10 2 cm. (16) 
Alternatively, for a given size of inhomogeneities a > a m ; n the scattering is strong for 

, . r~^ x f 30 cm for a = a mia , . 

A < A max = VaDd rs < ~ „ (17) 

|_ 10° cm lor a = a max 

A very important point is that this is a limitation of the wavelength from above: scattering is stronger 
for shorter wavelengths. This is in sharp contrast to the scattering in the interstellar medium, where the 
strength of the scattering increases at low frequencies. 

For the typical wavelength of observations A = 30 cm, both refractive and diffractive scattering occurs 
in a strong scattering regime with the total phase change (Eq. [Toj ) : 

( 1 diffractive ,_„, 
A H 100 refractive (18) 



Diffractive scattering is, though, only mildly strong. In comparison, in the ISM the small scale scintillations 
are saturated, while large scale are not. 
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The diffr active scattering angle 9]j will determine the lateral size a of the scatter-broadened spot: 

a « 6 D R LC « 5 x 10 7 cm (19) 

Interestingly, this is very similar to the size of the emission region of the Vela pulsar found using interstellar 
scintillations (Gwinn et al.1997) and consistent with other measurements of the emission sizes (Smirnova et 
al. 1996). 

Propagation of radiation through scattering media will also introduce temporal smearing of the pulse 
and related frequency decorrelation. Usually, these effects are due to the combination of three independent 
contributions ( Lee and Jokipii 1975b| ): (i) dispersive decorrelation effects arising from the nonzero d 2 ui/dk 2 
(which for sufficiently narrow receiver bandwidth effectively lead to a decrease in decorrelation bandwidth) 
(ii) diffractive decorrelation due to small scale scattering and (iii) refractive decorrelation due to the 
different transit times along the different ray paths. An important difference in our case is that dispersive 
decorrelation effects will be absent, since the waves are nondispersive. 

The extra path length introduced by strong diffractive scattering is 1 / '26 2 scat RLc ■ This will contribute 
to time delay and pulse broadening of the order 

t d = — — - =w 10 ^ sec, (20) 
2c 

(with corresponding decorrelation bandwidth Av£> — 2x 10 3 Hz). There is also a comparable group delay 
due to changing group velocity 

At group ~ 5— ps 10~ 5 sec. (21) 

In addition, the motion of the scatter-broadened spot due to both the rotation of the pulsar and the 
motion of scattering material inside the pulsar magnetosphere will produce refractive-type focusing and 
defocusing. The size of the diffraction pattern at Earth will be b ~ X/ddiff ~ 10 15 cm, where 

9 dlff ~0 D ^-^lO- 13 (22) 

is the visible size of the scattering spot (L is the distance of a pulsar to the Earth). The motion of this 
diffraction pattern due to the rotation of the pulsar will produce variations on a diffractive scintillations 
decorrelation time scale of 

r dl// -^P«^«10- 8 sec, (23) 
with associated decorrelation bandwidth 

Av diff = —^ «10 7 Hz. (24) 

27TT di // 

Motion of the scattering material inside the pulsar magnetosphere will produce diffractive scintillations on 
the same time scale. 

In addition the motion of the scattering disk of the size a will produce refractive type variations with a 
typical time scale 

tr ~ 9 R P = 10~ 4 sec (25) 

and refractive decorrelation bandwidth 

Av R fa 2 x 10 3 Hz (26) 
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2.2. Motion of plasma 

Plasma in the pulsar magnetosphere is moving along the magnetic field lines with relativistic velocities. 
This motion will affect the observed intensities in two ways: (i) motion of the inhomogeneities will produce 
temporal variations in the observed flux, (ii) we may expect to see aberration effects like relativistic 
contraction of the density inhomogeneities along the direction of motion which will produce a diffractive 
picture that is not circularly. The observed scattered image of the point source will have an elongated 
form with a ratio of axis w 7_l (7_i_ is a Lorentz factor of the motion in the plane of the sky - the parallel 
component of the perturbation's wave vector does not affect the scattering. For matter moving at an angle 
of about 45 degrees to the line of sight, we can estimate j± w y^y = 10. Thus, the image should be strongly 
elongated in the direction perpendicular to the direction of motion of matter. This effect, if resolved, will 
help to determine an absolute position of the rotation axis of the pulsar on the sky. 

To see how abberations appears in the general approach we note that for the case of frozen 
inhomogeneities, both motion of the plasma and relativistic contruction can be treated in the same 
systematic approach if we transform the distribution function of the index fluctuations from the plasma 
rest frame PkS(ui) to the laboratory frame. 

The Lorentz transformation of the power spectrum of density inhomogeneities P l , is related to the 
untransformed power spectrum P as 

P l (k l h k ± ) = P( 1 k h k ± ). (27) 

while the transformed correlation function of the power spectrum ip is related to the untransformed 
correlation function ip as 

iP l (r l h r ± )=^,r ± ). (28) 

(here the sign || denotes componets along the velocity of the medium). Using these relations it is easy to 
see that the relativistically moving medium is similar to a medium with anisotropic turbulence. 



2.3. Power-law distribution of inhomogeneities 

In turbulent media, it is common for the spectrum of the density perturbation to be a power law. 
Below we give order-of-magnitude estimates for this case (assuming that power law index a is less than 
four). 

There are some substantial differences between scattering by a screen with one scale and a power law. 
First, the correlation function of intensity fluctuations for the one-scale screen has scales of the order of the 
size of the inhomogeneities of the phase screen; while for the power law, the correlation function of intensity 
fluctuations has scales of the order of the Fresnel zone ^/XRlc- Secondly, the one-size screen can give 
strong focusing in the limit of scintillations (e.g., caustics can form if no diffractive size inhomogeneities are 
present); while in the case of the power-law screen, scintillation only approaches unity. 

The power law of density perturbations may extend from smallest scales a m i n = 1/qi ~ 5 x 10 2 cm to 
largest scales a max = 1/qo ~ 5 x 10 7 cm. To estimate the scattering for the power law spectrum of density 
perturbations ^N e oc k~ a , we note that the typical mean square density fluctuation due to waves with 
k w 1/a is SNg ~ N^(qo/k) a ~ 3 . The total phase delay due to the fluctuations on the scale a w 1/k for 
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strong turbulence SN^/N^ « 1 is then 

ID 



The mean scattering angle 



Atf, = JtLy/srtsl « D^qfr-WaT^X-H. (29) 
V a A 



= A^ = J^V^JwDVa^-^aV (30) 



is independent of frequency. For a < 4 (which includes the Kolmogorov spectrum a = 11/3), the scattering 
is mostly due to the small scale density fluctuations; while for a > 4, scattering is mostly due to the 
large scale refractive density fluctuations. The mean scattering angle (Eq.^) generally has only a weak 
dependence on the size (a- 1 ' 6 for the Kolmo gorov spectrum) . This suggests that the contribution from all 
scales may not be negligible. 

Estimates of the strength of the scintillations can be made for a power law distribution. The minimum 
scale that contributes to the strong scattering at a given frequency is (assuming a > 2) 

i 

a strong = — (-=:) " ' X~ 2 cm. (31) 



qo \D 

Moreover, for a given size a the scattering is strong for 



D S=2 5 

A<Amax=4/ — (q a) 2 w 20a 6 cm, (32) 

qo 



where we assumed go « 1/D and a = 11/3. 



2.4. Finite size of the emission region 

Untill now we neglected the finite size of the source in assuming that the initial unscattered wave 
corresponds to a plane wave. Finite size of the emission region may quench the scintillations. 

The criterion for scintillation to be quenched depends on whether the scintillation is weak or strong 
and whether the scintillation is refractive or diffractive. Weak scintillation, or strong refractive scintillation, 
will take place if the angular size of the source seen from the observer is less than the angular size of the 
scattering disk. For scattering very near the source (as magnetospheric scattering for a pulsar), scattering 
will take place if the linear size of the source is less than the linear size of the scattering disk. Strong 
diffractive scintillation will take place if the linear size of the source is smaller than the linear resolution 
of the scattering disk, seen as a lens. This can be expressed by requiring that the size of the source be 
less than X/8s, where 8s is the angular size of the scattering disk, seen from the source. For scattering at 
the light cylinder right near the source, 9s is the same as the total scattering angle 6d- Note that even if 
diffractive scintillation is quenched by finite size ofthe source, diffractive scattering will still lead to angular 
broadening. 

For a 0.1-sec pulsar, the light cylinder radius is Rl = 5 x 10 8 cm, the emission region size of the order 
of the polar cap radius Rpc ~ 1 km refractive scintillation will take place if the scattering angle is greater 
than AW « 2Rpc/Rl ~ 5 x 10~ 4 rad. Diffractive scintillation will take place if the scattering angle is less 
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than X/(2Rpc)- For 30 cm-wavelength observations of a 1-km polar cap, the required scattering angle of 
3 x 1CT 4 rad. 

For a msec pulsar, with the similar emission region size and cap and the light cylinder radius is 
Rl = 3 x 10 7 cm (P =5-msec) the refractive scintillation will take place if the scattering angle is greater 
than 1CP 2 rad while the limits on diffractive scattering will be the same as for 0.1-sec pulsar. 

Quenching by finite source size may also serve as an additional test to distinguish interstellar and 
magnetospheric diffraction: the latter should be more prominent in pulsars with narrow pulses. 



2.5. Millisecond Pulsars 

In this work we limited ourselves to fast rotating pulsars, for which the refractive index for radio 
wave propagation is independent of frequency. In slower rotating pulsars, the magnetic field near the light 
cylinder (adjusted for the relativistic motion of plasma) falls below the wave frequency; the plasma modes 
then have a "normal" dispersion (refractive index cx 1 — LUp/u 2 ). Scattering effects by such a plasma may 
change the simple frequency dependence (or independence) of the predicted effects. 

We expect that our order-of-magnitude estimates of scattering effects will hold for millisecond pulsars 
as well. Near the light cylinder, the parameter <5 (Eq. ^) is of comparable order for both millisecond 
pulsars and normal pulsars. The width of the scattering screen, which in millisecond pulsars occupies 
the whole magnetosphere, is also of the same order « 5 x 10 7 cm. The Fresnel radius is slightly smaller: 
ft = V ' Rlc^ ~ 4 x 10 cm. Quenching of diffractive scintillation may be more of a problem for millisecond 
pulsars, however, since they have larger emission beams. 



3. General considerations 

In this section we give general formulae for the scattering of electromagnetic waves by turbulent media 
within the pulsar magnetosphere. If the smallest size of the inhomogeneities is much larger than the 
wavelength, we can use the scalar wave equation to describe the wave scattering. Using the size of the first 
Fresnel zone rp = V^-Rlc ~ 10 5 cm for comparison purposes, we see that scattering by the inhomogeneities 
larger than rp occur in the geometric optics limit, while scattering by the inhomogeneities smaller than rp 
occur in the Fraunhofer limit. To avoid considering two models of scattering (as did Blandford & Narayan) 
we employ the Markov approximation, which is valid for short wavelengths satisfying both i) A a m j n 
(omin is the smallest scale of inhomogeneities) and ii) A <C l c where l c is the radius of coherence of the field 
(see below) (Rytov et al.). 

In the Markov approximation, the medium is assumed to be a collection of uncorrelated scattering 
phase screens - the correlation function of the refractive index fluctuations ip(r) is assumed to be delta 
correlated along the line of propagation: 

<A„(r) =< u(rx, z)u*(r' x , z') > = A(r ± - r' ± )8(z - z') (33) 

(Rytov et al.). Then A(r±) - the two dimensional correlation function - is related to the refractive index 
fluctuations as ^ 

A(rj_)= [ dzip n (r ± ,z)= f J (fc±r ± )P n (k ±) 0)4%. (34) 
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The physical meaning of the correlation function is that z = gives a length at which most of the energy 
of the wave is converted from ordered component into fluctuating. Related to the correlation function is the 
two dimensional structure function of the refractive index fluctuations for a homogeneous medium: 

D(r ± ) - 2 (A(0) - A(t ± )) = - f (1 - Mk ± r ± )) P n (k ± ,0)k ± dk ± . (35) 

Jo 

The scattering properties of the pulsar magnctosphcre are well approximated by that of an equivalent 
screen. The equivalent screen structure function of the phase fluctuations for a homogeneous medium is 



D 8n 2 r 2 r 4 r D-f 2 r qo 
D s (r±) = ^ D (^) = J dk ± k ± (! - Mr±k±))*AN e (k P er P ), (36) 



where r L = c/uo B is the Larmor radius. 

The normalized correlation function of the wave field u(r±, z) immediately after the screen is 

r(r ± , z) = exp{-L»(r ± ) s } =< u(rj_, z)u*(t ± , z) >= exp{--^D(rj_)}. (37) 

One can define several scales associated with the phase structure function D(r±)s: (i) the typical scale 
of D(p c )s = 1 determines the transverse coherence radius p c , so that the fields separated laterally by less 
that p c are strongly correlated; (ii) a longitudinal length z c = D ^ r± ^ which determines the minimal distance 
at which the Markov approximation is applicable; (iii) second (and larger) derivatives at zero displacement 
determine the typical small scale (rj_ -c a m in) behavior of the phase structure function 

2 _ d 2 D(0) s 

For larger scale inhomogeneities (such that a m j n ^> A), the transverse correlation radius after the 
screen is equal to the transverse correlation radius on the screen l± = p c ; while the longitudinal correlation 
length is lu = p 2 /\ > p c . 



3.1. Power law distribution 

We consider the power spectrum of the density perturbations (k) as the initial quantity for our 
approach. It is related to the the power spectrum of the refractive index perturbations P„(k) as 

P»(k) = * 2 %^. (39) 

We will use the following approximation for the refractive index perturbations: 

i?cxp{-fc 2 / gl 2 } 

rn[k) (i + k 2 /q 2 r/ 2 ' (4Uj 

where the outer scale is 1/qo = a max ss D, and the inner scale is l/qi = a m j n . The dependence of the 
various effects on the power law index is complicated. In what follows we assume that 3 < a < 4 (which 
includes the Kolmogorov value of a = 11/3). 
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The correlation function of the refractive index perturbations is 

The normalization of P„ is determined by the the mean square of the refractive index perturbations: 

(e 2 ) = Jd 3 kP n (42) 

8^r(f) /(^)\ 

s ^ ^xt^i)' l^rj (43) 

To normalize the spectrum of the refractive index perturbations, we assume that the turbulence is strong, 
i.e., the mean square of the relative density perturbations is of the order of unity: 



(«) 



Then, we obtain the following expression for the coefficient B: 

B = 87r3/2r (f) ,2 „ x 3 1Q 18 cm 3^ (45) 

?o 3 r(f-l) 

For a = 11/3, g = lO^cm" 1 , and <5 = 3 x 1CT 4 , 

Note, that for the chosen range of a, the normalization of the turbulent spectrum depends only on the 
outer scale; in this case, the increase of the power spectrum at small wave vectors outweighs the large phase 
volume at large wave vectors. 

For the power spectrum (Eq. [Io]) for distances much larger than the inner scale r± ^> 1/qi, the 2-D 
correlation function is 

A{V±) - 2^r( M +l) ' (46) 

where fi — a/2 — 1. 

For distances smaller than the outer scale r± <C qo and for a < 4 this simplifies to 
Bql 



A(rx) = 



27r(a - 2) 

Using the definition of B (Eq. [45| ) we find 

^ riV _Vy r(f) i 

A{r±) - qo 6 r(f-|)a-2 



(i - (W 2 ) 1 ^^ « and a < 4 (47) 



r(2-q/2) gap 2 
r(a/2) 1 2 j 



10 cm. (48) 



On the other hand, for scales much smaller than the inner scale r± <C 1/qi 

B 1o f-, 1 °--2A- a Jlt„. ow 4 " a 



^i) ri «i/ ?1 = 27r(a _ 2) ^ - 4«o «i - 2 ) r (^~) ) ^ « Vfc- (49) 
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So that the typical length scale za associated with A is 

A 2 



ZA 



A(Q) 



100cm. 



The 2-D structure functions of the index fluctuations are 

*J{rx) — n ( a Z 2 ) r(a/2) \ 2 ) *-/qi<^r± 
D(rx) = ^n^K' 2 ^ a r\ 
Using the definition of B (Eq. we find for l/q\ <C r± 

_ 8V^F 8 2 r(2-a/2) (qop\ a - 2 



rj_ < l/<7i 



q a- 2 T(f -§) V 2 7 



(50) 



(51) 



The condition D(r±)s = 1 defines /9 C , the correlation length: 



Pc = 



<7o 



8tt r(a/2)(a-2) 



n l/(a-2) 



BDP 9o 2 T(2 - a/2) 



2 



<7o a 



2T(f 



3^ 

2 2 / 



Dfc 2 0F <5 2 r(2-f) 



10 4 cm. 



(52) 



Note that in our case the correlation length is smaller than the Fresnel radius p c < ry. 

The other characteristic scale of D$, which determines the small scale r± <C l/qi behavior, is 



1 



dD s (r ± ) 



dr 2 



r ± =0 



BqfT(^) 



1/2 / x a/2 

(Zo 



10 b cm. 



(53) 



Thus the structure function is almost constant for r j_ <C 1 j q\ and has a typical correlation length p c given 
by Eq. ©. 

As the waves propagate through a turbulent medium, the angular spectrum will tend to become 
Gaussian. This occurs at (Lee & Jokipi 1975a) 



ZG 



(qi\ a r(3-a/2) 
Bq\ \q J r(2-a/2) 



10 9 cm. 



(54) 



This is much larger than the radius of the light cylinder, so that the angular spectrum does not become 
Gaussian in our considerations. 

There are several typical scattering angles associated with the power law distribution. First, there is 
the mean square angle: 



1 



-DV\D{ rA _) k=o| = 



47T 



r(- 



K ~ 2 lt a ~6xl0- 



(55) 



This is much smaller than the diffractive angle given in Eq. ( |14| ) for small scale inhomogeneities and is 
comparable to the refractive scattering angle Eq. (|l5|) obtained using simpler considerations. As we have 
mentioned, this is expected because our order-of-magnitude estimates for inhomogeneities dominated by 
small scales assumed that all of the power is concentrated at the diffractive scale. For the case of the power 
law, the power in the density perturbations is attenuated by the large phase space of long wavelengths 
refractive-type perturbations. 
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The other characteristic angle 9 C , which determines the typical width of the angular power spectrum 
for a power law distribution, is defined by the correlation radius 8 C = ^i-: 



2k 



BDk 2 q 2 r(2-a/2) 



8tt r(a/2)(a-2) 
The ratio of the two angles is 



l/(a-2) 



qo_ 

2k 



^Dk 2 5 2 T(2 - a/2) 
«b(a-2) T(^) 



l/(a-2) 



10" 



if 1/gi < p c \j%. 



0.2 for a 



(56) 
(57) 



The two angles scale differently with frequency: (0 2X ) is independent of frequency, while 9 C cx u<*-* 
(oc i^ 1 ^ 5 for Kolmogorov spectrum). They are just different measures of the angular spectrum; in principle, 
the two angles can serve to distinguish between Gaussian and power-law spectra of the density perturbations. 



3.2. Temporal broadening and Decorrelation 



In the absence of dispersive decorrelation, the intensity of the initial delta-pulse is a convolution of 
refractive and diffractive contributions only: 



P(z,t) = P D (z,t) * P R (z,t). 



The refractive Green function Pn(z,t) has a symmetrical Gausian form (Lee and Jokipii 1975b): 



PR(z,t) 



2ck 2 ( -2c 2 k 4 (t- z/vg) 2 



v / 2nA(0)D 



exp 



A(0)D 



(58) 



(59) 



with characteristic refractive decorrelation time and bandwidth 



TR = - 

c 



A(0)D 



1/2 



= 10" b sec 



27TTR 



2 x 10 4 Hz 



(60) 



The diffractive Green function Pr>{z 1 t) has an approximately exponential shape with a characteristic time 
of 

DO 2 

t d = w 10" 7 sec, (61) 



and an associated decorrelation bandwidth of 



2c 
1 

2txt d 



2 x 10 b Hz. 



(62) 



In case of power law spectrum the size of the diffraction pattern at Earth will be b f=a X/9diff ~ 10 17 cm, 
where 



10" 



(63) 



is the visible size of the scattering spot. The motion of this diffraction pattern due to the rotation of the 
pulsar will produce variations on a diffractive scintillations decorrelation time scale of 

b „ A 



T diff = J P 



C0i 



10 -10 sec 



(64) 
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with associated decorrelation bandwidth 

kvdi / / = — « 10 9 Hz. (65) 

2lTT diff 

These is a very larger decorrelation bandwidth. It is of the order of the observed frequency. Generally, the 
typical diffractive effects are weaker in the case of the power law spectrum of density inhomogeneities. This 
is due to the fact that in case of power law spectrum there is smaller power on the diffractive scales (given 
the fixed total power) than in the case of single size inhomogeneities. This makes the diffractive effects due 
to power law weaker and less likely to be observed. 

There though a possibility to distinguish power law from other types of density spectra. In case 
of power law distribution of inhomogeneities both to and 5vo are weakly dependent on the frequency 
(td oc i/ 2 / 5 ) because of the frequency dependence of 9 C . In contrast, if we have a Gaussian spectrum of 
density fluctuations, tr would be proportional to < 6 > and independent of frequency. This again may 
possibly serve as a test for the type of power spectra: for a Gaussian spectrum, the diffractive decorrelation 
time to should be independent of frequency; while for a power law spectrum, to oc zy 2 ( 4_Q )/( 2 ~ Q ). 

3.3. Correlations of intensities 

The correlation function of the wave intensities does not have an analytical representation (Rytov et 
al.). Asymptotic expressions are possible in the cases of weak or very strong scintillations. The division 
between these two limits is based on whether D{rj)s - the value of the structure function of the phase 
correlations at the first Fresnel zone - is <C 1 or 3> 1. This may be cast in the form of the ratio of the first 
Fresnel zone to the coherence radius V \D/p c : scintillations are weak for VXD/ p c <C 1, and strong in the 
opposite case. 

It is easy to see using Eqns ( p0| ) that in our case D(rf)s ~ 100 1. In this limit the scintillation 
becomes saturated with the scintillation index reaching 1. According to Prokhorov et al. 1975, there are two 

regimes for the spectrum of scintillations, separated by q* w l/V XDD(rf) « l/rj (this is equivalent 

to the separation between diffractive and refractive scintillation). For q <C q* ~ 1/r/ the modulations 
are weak, and the modulation spectrum is oc q 4 (Prokhorov et al. 1975 (4.40)) - these are diffractive 
scintillations. For q q* - refractive scintillations - modulations are strong and saturated: m z sa 1. 

If D(rj)s > 1 the scintillation index is (Prokhorov et al. 1975) 

m 2 «l+C%) s ^ (66) 

where C is of the order of unity. 

Thus we conclude that strong saturated refractive scintillations are possible inside pulsar magnetosphere. 

4. Observational Tests 

Here we summarize the predicted characteristics of the scattering inside the pulsar magnetosphere: 



diffractive scattering angle 
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diffractive scattering time 
diffractive decorrelation time 
refractive scattering angle 
refractive decorrelation time 
arrival time variations 



1(T 4 
1CT 8 
10~ 3 
10- 4 
1(T 4 



sec 



sec 



sec 



set- 



All these quantities are independent of frequency, but the strength of scattering increases with frequency. 
This frequency dependence of the strength of scattering is opposite to the interstellar propagation effects 
(which are weakest at larger frequencies). Thus the weak magnetospheric effects should be most prominent 
at high frequencies in nearby pulsars. 

In conclusion we discuss some observational facts that may possibly be attributed to magnetospheric 
scattering and propose future or follow up experiments to test the theory. 

Some of the propagation effects have possibly been observed. The most interesting and decisive (in 
our opinion) observations of Sallmcn et al. (1999) became known to us when most of the present work was 
completed. Two results of this work strongly support our theory: the frequency independent jitter in the 
arrival time, of the order 100/isec; and the frequency independent spread of the multiple components of 
the giant pulses, of tens of fisec, with various scattering times also of the order of tens of [isec. Using our 
results, these observations imply that the multiple structure of pulses is due to the multipath propagation 
inside the pulsar magnetosphere: it is frequency independent, but different rays propagate different lengths, 
aquiring a range of scattering times. 

Other observations that can be interpreted in favor of magnetospheric propagation include the 
following. Large sizes of the emission region(Gwinn et al. 1997, Smirnova et al. 1996 and Cordes et al. 
1983) may be due to diffractive scattering. Increase in the temporal broadening time td « 10~ 8 sec for 
nearby pulsars with dispersion measure DM < 20, relative to that extrapolated from more distant pulsars 
(Britton et al. 1998) may be due to diffractive decorrelation inside the pulsar magnetosphere. This td is 
comparable to our prediction (Eq. ^) if this is due to diffractive decorrelation inside the magnetosphere. 
Similarly, the predicted frequency independence of the decorrelation bandwidth (Eq. p4| ) naturally explains 
the results of Gwinn at al. (1999) - the diffractive broadening time td « 10 _8 sec will produce the observed 
5 MHz bandwidth gives. In the same manner the unusual decorrelation bandwidth of the PSR 0950 (Kramer 
et al. 1999) of 8 MHz is close to the predicted due to the diffractive scattering. If a follow up observation 
at a different frequency would show the same decorrelation bandwidth it will a strong argument in favor of 
the magnetospheric scattering. The enhanced intensity fluctuations (if compared with extrapolations from 
lower frequencies) at very high frequencies (30 GHz) found by Kramer et al. (1997) in several low DM 
pulsars maybe due to the scintillations inside the pulsar magnetosphere. 

The observations of B0950+08 at very low frequencies of 60 and 102 MHz by Smirnova & Shabanova 
(1992) which showed frequency independent narrowband variation of the pulse profile at both frequencies, 
with a characteristic bandwidth of 30 to 40 kHz may be explained as been due to refractive-type events 
with a typical refractive time tr~ 3 x 10~ 4 sec (Eq. p6| ) 

Experiments to detect effects of wave propagation inside the pulsar magnetosphere should use high 
frequencies for observation, and concentrate on nearby pulsars with low dispersion measure. Possible 
experiments will include a search for nondispersive (frequency independent) effects such as a time delay (as 
large as tens of microseconds) in the pulse arrivals, a diffractive decorrelation bandwidth of the order of 10 
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MHz, and microstructure periodicities (of the order of tens of microseconds) due to refractive scattering. 
Nearby strong pulsar, like PSR 0950, are best candidate for searches for magnetospheric effects. 

Other possible experiment that may be used to search for diffractive scattering inside the pulsar 
exploits interstellar scattering effects. Consider a thin interstellar screen which scatters pulsar radiation. 
The size of a patch on the screen is, approximately, the width of the pulsar diagram AW times the distance 
from the pulsar to the screen z. If there were considerable diffractive scattering inside the magnetosphere, 
the size of the coherent patch would be Xz/(0 s Rlc) <^ zAW. A smaller size of the coherent patch would 
change the observed properties of the diffraction pattern observed at Earth. Also, future high temporal 
resolution studies of pulsars (at time scales of nanoseconds) should provide more information on the weakest 
scattering events and possibly distinguish between interstellar and magnetospheric events. 

We would like to thank Carl Gwinn for numerous discussions and his valuable comments. AP thanks 
CITA for its hospitality during his stay in Toronto. This research was partially supported by NSERC. 
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A. Dispersion relation for streaming pair plasma 

If the average velocities of the electron and positrons of the secondary plasma are the same, the normal 
modes of strongly magnetized electron-positron plasma consist of three wave branches: extraordinary (X) 
and two coupled ordinary (O) and Alfven branches. Alfven mode cannot leave magnetosphere, while X and 
O mode may leave magnetosphere. 

For the forward propagating waves in the pulsar frame we have 



7p ^(l-, p co s9 /c)^ _ / fee (l - j^r) , if « £ 



u x =kc(l- W^ose/c) ) 



V U B / IP 



i ( 1 7 P ai (t> p /c-cos0) u sin 8 \ I , / 4 7 w 2 s in 2 8 \ t / 2 

= kc[l *— ^ + ^ c2k2 { l_ VpCOB9/cy2 ) = < kc (1 + j if uj/j^b « 9 « l/ 7p if 



/ ,2 2 2 2 \ f kc cos 6> fl - , " p 2 ) if < " p 

c A = fcccosfl ( 1 - T ^ - C J sin .f ) = I ) *f?' 2 N ^ if 
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